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Abstract — A post-processing adjustment technique 
which aims for enhancement of dual-frequency SURF 
(Second order UltRasound Field) reverberation-noise 
suppression imaging in medical ultrasound is analyzed. 
Two variant methods are investigated through 
numerical simulations. They both solely involve 
post-processing of the propagated high-frequency 
(HF) imaging wave fields, which in real-time imaging 
corresponds to post-processing of the beamformed 
receive radio-frequency signals. Hence the transmit 
pulse complexes are the same as for the previously 
published SURF reverberation-suppression imaging 
method. The adjustment technique is tested on 
simulated data from propagation of SURF pulse 
complexes consisting of a 3.5 MHz HF imaging pulse 
added to a 0.5 MHz low-frequency sound-speed 
manipulation pulse. Imaging transmit beams are 
constructed with and without adjustment. The post- 
processing involves filtering, e.g. by a time-shift, in 
order to equalize the two SURF HF pulses at a chosen 
depth. This depth is typically chosen to coincide with 
the depth where the first scattering or reflection 
occurs for the reverberation noise one intends to 
suppress. The beams realized with post-processing 
show energy decrease at the chosen depth, especially 
for shallow depths where in a medical imaging 
situation often a body-wall is located. This indicates 
that the post-processing may further enhance the 
reverberation-suppression abilities of SURF imaging. 
Moreover, it is shown that the methods might be 
utilized to reduce the accumulated near-field energy 
' of the SURF transmit-beam relative to its imaging 
region energy. The adjustments presented may 
therefore potentially be utilized to attain a slightly 
better general suppression of multiple scattering and 
multiple reflection noise compared to for non-adjusted 
, SURF reverberation-suppression imaging. 
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I. Introduction 

THE SURF imaging synthetic transmit-beam for 
reverberation suppression, below referred to as the 
SURF beam, is generated from the difference between 
two subsequent ultrasound imaging pulses transmit in the 
same direction [T] . Each of the transmissions consists of 
a dual-frequency pulse complex with a conventional high- 
frequency (HF) imaging pulse added to a low-frequency 
(LF) tissue manipulation pulse. The high frequency /h 
may typically be around 10 times the low frequency /l- 
The LF polarity is switched for the second transmis- 
sion so that the LF pressure experienced by two HF 
pulses becomes of opposite polarity. Therefore, due to the 
nonlinear effect of pressure-dependent compressibility of 
tissue, the LF-pulse dependent effective speed of sound 
experienced during propagation of these two HF pulses 
differs. As the covered distance increases, a difference is 
thus accumulated between the propagation time required 
for the HF pulses of the two transmissions. The difference 
obtained when subtracting the two propagated HF pulse 
components hence grows with covered distance until they 
possibly are separated by a delay corresponding to a 
180° phase-shift. This effect may be utilized during image 
reconstruction in order to suppress multiple scattering 
and multiple reflection artifacts (in the rest of this paper 
considered as and denoted as reverberation artifacts), as 
described and analyzed in Ref. (Ij where also combined 
SURF and pulse-inversion image reconstruction is intro- 
duced. In the citet paper, SURF transmit beams are 
simulated for two different pulse and transducer setups 
and compared to conventional fundamental and pulse- 
inversion second-harmonic transmit-beams. In this paper, 
the theory behind SURF reverberation suppression is 
further described in Section [TTl For a recent review of tissue 
nonlinearity, second-harmonic tissue imaging, and related 
research, see Ref. i. 

The reverberation suppression ability of SURF imaging 
is illustrated by its synthetic transmit beam being re- 
duced near the transducer, as is also the second-harmonic 
transmit-beam in tissue harmonic imaging ^ . 

Dual- frequency pulses have also been used for ultra- 
sound contrast-agent microbubble detection . Within 
this context, SURF imaging is sometimes denoted radial 
modulation imaging. The contrast agent detection has 
been shown in vivo [1^. In Ref. 0, the versatility of 
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SURF imaging is illustrated, and fundamental concepts 
in relation to the use of dual-frequency pulse complexes 
in ultrasound imaging are investigated both regarding 
nonlinear propagation and nonlinear scattering. A variant 
of the time-delay post-processing adjustment analyzed 
in the present article is also considered. Furthermore, 
dual-frequency pulse insonification has been employed for 
imaging the nonlinearity parameter of different media [12- 

The purpose of the adjustment methods analyzed in the 
present paper is to further improve SURF reverberation 
artifact suppression, especially if the dominating rever- 
beration noise is due to a strong first scattering/reflection 
taking place around a depth Za in the near-field. Such 
suppression would be helpful in clinical applications e.g. 
where some weakly scattering anatomic feature is to be 
studied within a hypoechoic region. Then strong rever- 
beration noise artifacts risk to camouflage the features 
intended to be visualized. 

The transmitted and propagated pulse complexes uti- 
lized are equal to those in non-adjusted SURF imaging. 
As described in Section |TT1 a transmit synthetic SURF 
HF wave field is generated as a difference, with one of 
the propagated HF wave fields being adjusted, before the 
adjusted SURF difference is calculated. Subsequently a 
transmit beam is generated in a conventional way from the 
the root-mean-square (RMS) or the temporal maximum of 
the adjusted difference field. 

The rest of this article is organized as follows: First 
the SURF reverberation-suppression imaging method is 
reviewed. Then comes a description of three different kinds 
of time-shifts of importance for conceptual understanding 
of the here introduced post-processing methods. Two post- 
processing signal adjustments are depicted: adjustment 
by a time-shift and adjustment by a general filter. The 
introduction of the general filter, which in addition to 
time-shifting also adjusts the pulse form, is motivated 
by the possibility that the two HF wave forms may be 
distorted if the experienced LF pressure varies along the 
HF pulse at a given spatial location. This variation is of 
opposite sign for switched LF polarity. 

Subsequently, a SURF beam gain-factor associated to 
the time-shift adjustment is explained. Then is defined a 
specific beam quality ratio Qz^ utilized to score transmit 
beams regarding suppression of reverberation noise stem- 
ming from a first scattering/reflection emerging at a depth 

The findings in the Results section are based on post- 
processing adjustment of data generated by nonhnear 
propagation simulation of two transmit SURF pulse com- 
plexes of opposite LF polarity. The transducer and pulse 
complex setup is described, then the effect on the spatial 
development of the gain factor due to the post-processing 
time-shift is shown. Then the beam-quality ratio of the 
calculated beams is provided. Finally, SURF beams are 
constructed and compared for non- adjusted data as well 
as data adjusted by time-shifts and by general filters. 



IL Theory and Methods 

A. SURF reverberation-suppression imaging outline 

This subsection describes SURF reverberation- 
suppression imaging without post-processing adjustment. 
Each SURF transmit pulse is a dual-frequency complex 
composed of a high-frequency (HF) imaging pulse added 
to a low-frequency (LF) tissue manipulation pulse. 
Reverberation noise suppression is achieved by taking 
advantage of an accumulated propagation time-delay of 
the HF pulse, introduced by the LF pulse. In the medical 
ultrasound regime, the nonlinear LF pressure- dependent 
speed of sound experienced by the HF pulse is 11 

c(r) =co(l+/3n/^PL(r)), (1) 

where cq is the speed of sound with no LF manipulation 
present, Pn is the material nonlinearity parameter, k the 
compressibility and Ph{f) is the LF pressure at the spatial 
location r. The effective speed of sound when propagat- 
ing in a positive LF pressure, and therefore compressed 
medium, thus becomes higher than in case of the medium 
being expanded due to a negative LF pressure. 

For SURF reverberation suppression image reconstruc- 
tion, one transmits two SURF pulses for each transmit 
focal zone and direction. The HF pulse is equal for 
both transmissions, while the polarity of the LF pulse 
is switched for the second pulse. These HF pulses thus 
propagate at different effective speed of sound due to 
the opposite polarity of the associated LF pulses. Conse- 
quently a relative propagation time-shift r(r) is generated 
between them, which is small at shallow depths, but is 
accumulated during forward propagation within the com- 
pressed/expanded medium. The synthetic SURF transmit 
beam is constructed from the difference between the HF 
wavefields. For a given r, its maximum amplitude is for 
the accumulated time-shift between the two HF pulses 
r(f) = 1/(2/h). 

Due to the amplitude reduction of all reflected or 
scattered parts of a wave, an imaging pulse experiences 
negligible sound speed manipulation after being scattered, 
since the scattered or reflected LF wave it propagates 
in conjunction with is significantly weakened. Therefore 
the multiply scattered or reflected HF wave accumulates 
a smaller relative time-shift than the direct forward- 
propagating part. The reverberated parts of the receive 
HF pulses thus give negligible contribution to the receive 
SURF difference signal if the first scattering/reflection 
takes place near the transducer where the accumulated 
time-shift r(r) during forward propagation is small [l[. 
The cited paper describes three different classes of re- 
verberation noise, and the simulated beams are used to 
compare the reverberation suppression abilities of the 
respective imaging methods. This is done by comparison 
of an introduced quality parameter Q (also stated here in 
Section HV|). which measures the ratio between the beam 
energy ratio within the imaging depth region and the beam 
energy closer to the transducer, thus quantifying how well 
Class I and II reverberation noise is suppressed. 
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B. Time-shift terminology 

There are three tirae-shifts of importance for for con- 
ceptual understanding of the presented post-processing 
adjustment method: (i) tq between the center of the HF 
and LF pulses in a transmit SURF pulse complex, see 
Ref. [l|,[7| for further elaboration, (ii) r(r) accumulated as 
the HF pulse propagates. It is due to the nonlinear sound 
speed change ([1]) induced by the LF part of the SURF 
pulse complex, as explained above in Section IH-AI (iii) 
the post-processing reverberation-suppression adjustment 
time-shift r^, applied to one of the propagated HF fields 
before forming the adjusted SURF difference field, as 
described below in Section Ill-Ci 

C. SURF transmit-beam adjustment by a post-processing 
time-shift Ta 

In the model provided below, s+(r, t) represents the 
HF part of a SURF wave field propagated with positive 
LF manipulation polarity, while s_(r, represents the 
corresponding field in case of negative LF polarity. In this 
model, these HF wavefield components differ solely by a 
spatially varying accumulated time-shift r(r) induced by 
the LF sound speed manipulation. In addition, a con- 
stant post-processing time-shift is applied to s_(r, t) 
to generate the field denoted s_(r, t). The HF field after 
propagation without LF manipulation is denoted s(r, t). 
Then s+(r,t) = s(r, t + r(r)/2) and s_(r,t) = s{f^t — 
T{r)/2-Ta). The time-shift adjusted SURF HF difference 
field then becomes: 

SA{r,t) =s^{r,t) - 5_(r,t) 

=s (r, t + r(r)/2) -s{r,t- r(r)/2 - r^) . (2) 

For a non-adjusted field on the s(r, t) = 9^e{5(r, t)e*^°*} 
form, where ^(r , t) is a complex pulse envelope and ujq the 
center angular frequency, the difference hence becomes 

5A(r, t) = 9lt{ + s(r, t + r(r) /2) e^^^^^+^^^^^Z^) 



s{f,t- r(r)/2 - Ta) e 



ioJo{t-r{r)/2-Ta) 



}■ (3) 



Having the envelope x unchanged during a time-shift, 
which is a narrow-band approximation, yields: 



SA 



Jujot \ iuJoT{f)/2 



-iu;o(r(r)/2+ra)l 



=lHe{s(r,t) 
= 2 sin[c^o(r(r) + ra)/2] 3m {s{t)e'^'' } , 



} 



(4) 



= G(cJo, r(f) +ra) 

where G may be interpreted as a spatially varying gain 
factor regulating the amplitude of sa(^, relative to the 
non-shifted field ^(r, t). The temporal maximum of sa may 
thus be estimated from knowledge of r(r), r^, and the 
temporal maximum of s(r, t) without LF manipulation: 

G{uJo,r{r)^ra)3m{s{r,t)e'^''} 
= \G{uJo,r{r) + r^)! max bm {s(r, t)e'^^'} 

or equivalently for the RMS. The first maximum and 
minimum of |G(co'o, r(r) + ra) \ are thus {2, 0} for uJo{r{f)-\- 
ra) = KO}. 



max 

t 



SA 



max 

t 



With no time-shift adjustment applied, G = for z = 
(at the transducer), while application of adjusts the 
depth of maximum suppression. The time-shift adjustment 
Ta giving maximum suppression of sa(^, ^) at the on-axis 
depth Za may be predicted from (|4]) for G{uJo^r{za) + 
Ta)) = 0, which equally gives s^{za,t) = S-{za,t). 

D. SURF transmit-beam adjustment by a general post- 
processing filter 

The nonlinear LF manipulation might generate HF 
waveform change differing between s+(r, and S-{f^t) 
in a manner that cannot fully be described by a time- 
shift, e.g. due to pulse compression or expansion. Then 
a time-shift adjustment Ta{za) does not exist to make 
5+(r, = s-{f,t) at z = Za. However, instead of a pure 
time-shift a more general filter f)^^ may be applied to s(r, t) 
to make 5+(r, t) = l)^^ {s_(r, t)}(r, at the reference 
depth Za. Then the SURF difference field becomes 



SA{r,t) = 5+(r,t) - l)^Js_(r,t)}(r,t). 



(5) 



In this work, Ijza is chosen as the Wiener filter making 
5+(r, t) = s_(r, t) on-axis at the reference depth Za. This 
filter is then applied to the whole s_(r, t) field. 

E. Specific beam quality ratio 

For reverberation suppression in practical imaging, it is 
important that the sensitivity to multiple scattered and 
multiple reflected waves is low compared to the sensitivity 
to direct back-scattered waves from the imaging depth. 
Thus reduction of the transmit beam at the depth Za by a 
factor C is of no use if the beam is equally suppressed by C 
within the imaging region. To estimate the reverberation 
suppression characteristics of SURF transmit beams in 
case of a first scattering/reflection at the depth Za, a 
specific beam quality ratio Qz^ is here introduced as the 
transmit beam energy within the imaging region divided 
by the summed beam energy at the depth Za' 



(6) 



z^Zr, r=0 6'=0 



r=0 9^0 



where E{f) is the beam energy at the spatial coordinate 
f described by the cylindrical coordinates {z^r^O). The 
imaging depth region is within z G [zn^Zf] and is equal 
to the focal region given the active transmit aperture and 
imaging frequency. 

III. Results 

A. Simulation setup 

The raw data utilized to generate the results presented 
below is the same computer simulated fields of HF pulses, 
propagated in conjunction with LF manipulation of pos- 
itive and negative polarity, as generated for the 3.5 MHz 
setup in Ref.UI. These HF wave fields are used as s + (r , 
and (non- adjusted) 5_(r, t). The two adjustment variants 
(the pure time-delay filter and the delay and pulse-form 
filter) are applied to 5_(r, t) to generate s_(r, t). The 
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Figure 1. Transmit aperture. Central region: HF and LF, outer 
region: only LF. 

modeled transmit aperture is axisymmetric consisting of a 
central region with both HF and LF transmission and an 
outer region where only LF is transmitted. Fig.[T]shows the 
geometry of such HF and LF transmission regions. The 

Table I 

Transmit pulse parameters for SURF beam generation. 



/h 


/5H,6dB 


/l 


/5L,6dB 


POH 


Pol 




[MHz] 


% 


[MHz] 


% 


[MPa] 


[mm] 


[mm] 


3.5 


50 


0.5 


25 


3.5 


0.85 


7.1 








TO 




^/ 






[MPa] 


mm] [mm 


[MS] 


[mm] 


[^s] 






10 


82 82 


-0.2 


60 


130 





/ho is the transmit HF frequency, (3 the fractional frequency band- 
width, ^0 transmit surface pressure, ao the outer aperture radius, F 
the focal depth, and /h the imaging frequency. H and L refer to high 
and low frequency, ro is the delay between HF and LF transmission. 
The imaging region is 2; G [zn, Zf]. 

HF and LF frequencies, bandwidths, excitation pressures, 
and aperture radii are described in Tab. HI where also the 
chosen delay tq between transmission of the HF and LF 
pulses are given. The apertures are modeled without kerf. 

B. Estimated time-shifts and gain factors 

Fig. [2] displays the estimated time-shift r(r) between the 
simulated HF fields 5+(r, t) and s_(r, t), as well as the 
corresponding gain factor estimates, calculated from the 
definition in Eq.lH The time-shift adjustment by corre- 
sponds to a vertical movement of the r{z) curve, which 
infiuences the gain factor magnitude, especially within 
the near-field. To illustrate this. Fig [2] in addition shows 
r{f) -^Ta and G(a;o,r(r) + r^) for Ta = {-21, 7.1, 54} ns, 
corresponding to {—27°, 9°, 69°} phase-shifts at 3.5 MHz. 

The depth-dependent optimum time-shift Ta{za)^ with 
Za as the chosen depth of increased suppression, is found 
by maximizing Qza{^a) under variation of r^. Fig. [3] dis- 
plays the resulting relation between this depth and the 
time-shift. 

C. Adjusted synthetic transmit imaging beams 

The s-{f^t) field is adjusted both using the optimum 
time-shift ra{za) and using the general filter l)^^. The 
SURF difference HF field is formed for both adjustment 




10 20 30 40 50 



Figure 3. Optimum time-shifts Ta as function of the suppression 
depth Za- 




5 10 15 20 25 30 35 
2 [mm] 



Figure 4. Specific beam quality Qza- '• 'Ta adjusted, : 

f)2^ adjusted, : non- adjusted, — — — : HF beam without LF 

manipulation (fundamental imaging). 

methods. The depth choice of optimum suppression is var- 
ied within Za G [1,55] mm, and resulting adjusted SURF 
transmit fields for Za = {5, 10, 20, 30, 40, 55} mm with 
5_(r, unmodified, modified by the optimum time-shift, 
and modified by l)^^, are shown in Figs. [BHTTl displaved at 
the end of the paper. 

D. Beam quality ratios 

Fig. m shows a comparison between the specific beam 
quality ratio Qzai^a) calculated for beams with different 
adjustments, for the non-adjusted beam, and for the con- 
ventional fundamental imaging HF transmit-beam with- 
out LF manipulation. 

IV. Discussion 

As illustrated in Figs. MM the realized time-shift ad- 
justments on the simulated data generate a decrease of 
on-axis amplitude of the individually normalized SURF 
beams of ~{10, 9, 4, 3} dB at the chosen reference depths 

= {5, 10, 20, 30} mm of increased suppression. 

The time-shift adjustment of 5_(r, t) equals a move- 
ment of the on-axis zero-crossing depth coordinate of the 
r{z) + Ta, and consequently a change of depth where the 
gain-factor model, given in Eq.|4] predicts G = and hence 
zero SURF beam amplitude. However, the realized time- 
shift adjusted SURF beams display a beam amplitude at 
Za which gets less suppressed with increasing Za. An ex- 
planation to this is that the effect of the LF manipulation 
is not solely a pure time-shift but also involves some pulse- 
form distortion, e.g. due to pulse compression or expan- 
sion. This distortion may be caused by an experienced LF 
pressure which varies along the time axis for a HF pulse 
having propagated to a given spatial coordinate. 
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(a) Ta = 0. Horizontally: (b) Ta = 7.1ns. Axes as in (a). 

depth, vertically: lateral 

coordinate. 




(c) Ta = 54ns. Axes as in (a), (d) On-axis r(z) + Ta and 

corresponding G(z). : 

Ta = -21 ns, — — — : Ta = ns, 

Ta = 54 ns. 

Figure 2. Top sub-panes: Estimated time-shift t(z) + Ta [ns] between s+(r, t) and s_(r, t). Bottom sub-panes: corresponding gain factor 
estimates G(z) [dB]. 



The shown adjustments using a general filter () generate 
full transmit beams suppression at the Za depths on-axis. 
However, the corresponding SURF difference HF pulse 
gets significantly distorted within the imaging region when 
20^^a^30mm, as exemplified for Za = 30 mm in Fig. [9l 
where the on-axis pulses are significantly less concentrated 
in time within the imaging region than for Za<20mm. 
This effect is due to the edge waves arriving on-axis in 
this region, hence enforcing l)^^ to be complicated in 
order to absolutely equalize s+(r, t) and s_(r, t). Longer 
difference pulses in the imaging region cause decreased 
range resolution. 

The determination of f);^^ is here based solely on the on- 
axis fields. Hence the adjusted transmit-beam suppression 
may be insufficient off-axis. Optionally the calculation of 
could be based not only on the on-axis, but also the 
off-axis characteristics of 5+(r, t) and 5_(r, t). However, 
the adjusted transmit beams realized here and partially 
shown in the results section prove to have little energy off 
the beam axis for adjustments with Za<35 mm. 

The increased off-axis energy with increasing Za is re- 
flected in Fig. S] where Qz^ is shown to decrease from 
~75dB at = to '-«^OdB at Za = 30 mm for l)^^ 
adjustment. The corresponding time-shift adjusted beams 
give a Qz^ decreasing from ~30dB to OdB within the 
same depth interval, while the non- adjusted SURF beam 
at any Za depth gives a Qz^^lOdB below what is attained 
using the time-shift adjustment. There is no significant 
gain in Qz^ by utilizing l)^^ adjustment instead of time- 
shift adjustment for ^^^20 mm. 

For the adjustments realized on the example dataset, the 



beams for Za^SOmm suffer from giving high beam energy 
within the near-field. Thus the reduced sensitivity to 
reverberation-generating scatterers at Za is then attained 
at the cost of high sensitivity to reverberation-generating 
scatterers within the near-field. Since an imaging situ- 
ation with strong reverberation-generating scatterers at 
great depths is likely to also have strong reverberation- 
generating scatterers within the near-field, the method 
thus has limited benefit compared to non- adjusted SURF 
when 2:^^30 mm for the example dataset. 

Determination of l)^^ in a real imaging situation might 
be facilitated if d priori information is known about how 
s+(r, and s_(r, t) differ depending on Za- Potential 
sources for gathering of such information are: (i) received 
directly back-scattered HF signals from the Za depth, (ii) 
pre-recorded hydrophone measurements of the transmit 
waves at Za^ and (iii) numerical transmit wave field simu- 
lations, as utilized in this work. Based on either of these 
information sources, a set of pre-defined l)^^ could be 
generated from within the scanner operator may choose 
at will. 

Fig. [2] shows that application of a time-shift adjustment 
of opposite sign than for movement of the maximum 
suppression depth Za into greater depths, the gain factor 
instead grows within the whole near-field and the whole 
imaging region. The signal magnitude within the focal 
region is thus increased at the cost of less reverberation 
suppression. This may be favorable for the general image 
quality in case of little reverberation noise, and may be 
utilized to improve the sensitivity and signal-to-noise ratio 
(SNR), or may alternatively be traded-off into increased 
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Figure 5. General beam quality ratio Q as function of Ta for the 
time-shift adjustment (top pane), and Za for the filter adjustment 
\)za (bottom pane). 

penetration depth. This is vahd for situations including 
the one presented here, where the accumulated time-delay 
within the imaging region is smaller than the gain factor 
maximizing delay 1/(2/h). 

As shown in Fig. [2l the relative change in gain-factor 
magnitude is larger within the near-field than within the 
focal region, for small time-shift adjustments r^. The 
signal model predicts this effect to be amplified due to the 
lesser relative sensitivity to changes in the argument of the 
sinusoidal which models the gain factor, Eq.Sl in case of 
G being close to its maximum than when it is small. This 
observation inspires to assume that Ta or l)^^ adjustment in 
SURF beam generation might be useful not only to better 
suppress reverberations arising from a single depth, but 
also to get improved general reverberation suppression by 
decreasing the ratio between the total energy within the 
near-field z G [0, Zn] and the energy within the imaging 
region z G [zn-,Zf]. This is illustrated in Fig. [5l which 
shows the respectively \)z^ infiuence on the general beam 
energy quality ratio Q, as defined in [1]: 

2/ CO 27r / Zn oo 27r 

z=zn r=0 6'=0 / 2=0 r=0 6'=0 

where E{r) is the beam energy. The attained ratio Q, as 
displayed in Fig. O has a maximum around the time-shift 
Ta ^ 12.5 ns or the general filter \)za^ib mm-, thus sug- 
gesting that such adjustments of S-{f^t) generate SURF 
beams generally more advantageous for reverberation sup- 
pression than without adjustment. The increase in general 
beam quality ratio is a modest '-«^1.5dB. 

V. Concluding remarks 

This work indicates that the SURF reverberation sup- 
pression transmit-beam may be adjusted to increased 
canceling of reverberation noise generated due to first 
scattering/refiection from a certain depth of choice Za- 
The adjustment is done without change of the transmit 
SURF pulse complexes utilized, but instead by application 
of some filter, e.g. a time-shift, on one of the propagated 
HF wave fields before SURF difference calculation. Dif- 
ferent filters/time-shifts are utilized for different depths 
of first multiple scattering/refiection of the reverberation 
noise desired to suppress. The analysis of the Class I and 
II reverberation noise suppression abilities are based on 
synthetic SURF transmit-beams after adjustment. 



If the depth position Za of this first scattering/refiection 
to be suppressed is outside of the near-field, the reduction 
of sensitivity to scatterers at Za is attained at the cost of 
increased sensitivity to scatterers within the near-field. 

A strength of the method is that, as long as the re- 
ceive HF imaging signals for propagation with both LF 
polarities are stored, there is no need for re-transmission 
of pulses to regulate the region of maximum reverberation 
suppression. The user may thus on-the-fiy regulate the 
adjustment ad hoc^ while interactively observing which 
adjustment that generates the most clear image for the 
imaging region of interest depending on the nature of 
the reverberation noise present, or optionally to increase 
the sensitivity at large depths. The choice of adjustment 
may also be regulated to improve saved images as long 
as the HF receive IQ- or radio-frequency-data needed to 
form the SURF difference image is stored. The method 
thus gives a larger fiexibility regarding regulation the 
reverberation suppression compared to e.g. pulse inversion 
imaging where normally a new pulse couple has to be 
emitted if one aims to modify the transmit beam. 

The realization of the methods presented is limited to 
computer simulations of pulses in a homogeneous medium. 
Additional simulations for an inhomogeneous medium, as 
well as experimental setups with a real transducer and 
phantom, and in vivo studies, are natural steps to further 
validate the feasibility in real imaging situations. 
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(a) SURF pressure beams. 
Left: HF difference with Ijz^ 
adjustment, middle: HF 
difference with Ta adjustment, 
right: HF pressure 
non-adjusted. 
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(c) On-axis SURF pulses. 
Left: HF difference with Ijz^ 
adjustment, middle: HF 
difference with Ta adjustment, 
right: HF pressure 
non-adjusted. 



(b) On-axis max. SURF HF 
difference pressure. Top: 
individual normalization, 
bottom: common norm. 
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(d) On-axis SURF pulses. Top 
pane: at the depth Za, 
bottom: at the focal depth. 



Figure 6. Transmit beam and on-axis pulse comparisons, 
: HF pressure non- adjusted, — 



Za = 5 mm 



for |(b)| and 1(d) 
adjustment Ij^ 



SURF pulses and fields generated from simulated data. Line notation 
HF difference with time-shift adjustment, : HF difference with filter 
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